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Abstract 

Karma and Rappel [1] recently developed a new sharp interface asymptotic analysis 
of the phase-field equations that is especially appropriate for modeling dendritic growth 
at low undercoolings. Their approach relieves a stringent restriction on the interface 
thickness that applies in the conventional asymptotic analysis, and has the added ad- 
vantage that interfacial kinetic effects can also be eliminated. However, their analysis 
focussed on the case of equal thermal conductivities in the solid and liquid phases; when 
applied to a standard phase-field model with unequal conductivities, anomalous terms 
arise in the limiting forms of the boundary conditions for the interfacial temperature 
that are not present in conventional sharp- interface solidification models, as discussed 
further by Almgren [2]. In this paper we apply their asymptotic methodology to a 
generalized phase-field model which is derived using a thermodynamically consistent 
approach that is based on independent entropy and internal energy gradient function- 
als that include double wells in both the entropy and internal energy densities. The 
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additional degrees of freedom associated with the generalized phase-field equations can 
be chosen to eliminate the anomalous terms that arise for unequal conductivities. 

PACS: 81.10.Aj, 64.70.Dv, 66.10.Cb, 81.10.Mx 
Keywords: phase-field model, sharp interface analysis, thin interface analysis, dendritic 
growth, isothermal asymptotics. 


1 Introduction 


Various asymptotic limits can be used to relate diffuse interface models of solidification 
to conventional free boundary problems that involve sharp interfaces of zero width (see, 
for example, [3]). In particular, the size of the diffuse interface width i relative to other 
characteristic length scales in the problem determines the specific form of the free boundary 
problem that is obtained as i becomes small. In this paper we will be especially concerned 
with the effects of capillarity and interface attachment kinetics, as described by the modified 
Gibbs-Thomson boundary condition, 


T-Tm 

Ly/cp 


-m - v/v 0 , 


(i) 


which relates the interface mean curvature and normal velocity V, to the dimensionless 
interface temperature, u = (T — T^) /[Ly /cp], where T is the temperature, Tm is the melting 
point, Ly is the latent heat per unit volume, and Cp is the heat capacity per unit volume. Here 
£ c — TmJ /[Ly/cp] is a capillary length based on the surface free energy 7, and Vo — n[Ly/cp] 
is a characteristic velocity based on the kinetic attachment coefficient, fi. For discussion 
purposes we consider the growth of a dendrite with tip radius R into an undercooled melt 
[4]; in more general circumstances, R might be representative of a typical macroscopic length 
scale such as a container size. For dendritic solidification at large undercoolings the growth 
is rapid and the radius of curvature of the tip is relatively small, so that the effects of 
capillarity and kinetics on the interface temperature in Eq. (1) can be significant. In this 
regime, sharp interface limits of the phase field equations have been performed [5-10] that 
assume the interface temperature u is of order unity in the small parameter £/R, which allows 
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the contributions from capillarity and kinetics to be of order unity as well. In this limit l is 
also assumed to be small compared to £ c , which presents a stringent resolution requirement 
for a numerical computation that aspires to describe this limiting case. At low undercoolings, 
on the other hand, dendrites grow more slowly and have a larger radius of curvature, so it 
is reasonable to model the effects of capillarity and kinetics as small corrections. Karma 
and Rappel [1] refer to the corresponding analysis as the “thin interface limit,” wherein one 
assumes Icli but allows £ ~ l c . Almgren [2] has described this analysis as “isothermal 
asymptotics,” since to leading order in l/R the temperature is isothermal throughout the 
interfacial region, with u — 0(£/R). A technical difficulty with this analysis is that for 
unequal thermal conductivities in the liquid and solid, the thin interface limit can lead to 
a discontinuity in the temperature field across the interface as well as the appearance of 
an anomalous term proportional to the interfacial temperature gradients in the modified 
Gibbs-Thomson equation. Almgren shows that a continuous temperature field and a Gibbs- 
Thomson equation independent of the interfacial temperature gradients can be achieved by 
appropriate alterations of the interpolation functions that are used to represent the bulk 
thermophysical quantities across the diffuse interface. 

The standard phase-field model, a version of the Model H of Halperin, Hohenberg, and 
Ma [11], was originally derived on the basis of a gradient functional for the free energy of the 
system, with dynamics that ensure that the free energy decreases with time [12-16]. The 
time rate of change of the order parameter is coupled to the energy equation to incorporate 
the release of latent heat as solidification occurs. Penrose and Fife [17] (see also [18]) derived 
the phase-field equations on the basis of positive entropy production for a system that is 
based on a gradient functional for the entropy of the system. Here we consider a generalized 
phase field model that is based on gradient functionals for both the entropy and internal 
energy of the system; such models have been considered by a number of authors [19-23]. 
The generalized model contains additional degrees of freedom that we use here to eliminate 
the anomalous terms that arise in the asymptotic analysis of Karma and Rappel in the case 
of unequal thermal conductivities. 
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The kinetic attachment coefficient /j that appears in the characteristic velocity Vo in 
Eq. (1) has been measured for relatively few materials. For this reason, and because the 
effect of kinetics is usually expected to be small at low to moderate rates of solidification 
for metals and other molecularly-simple materials, many sharp interface models ignore the 
effect of kinetics by setting l/fi, — 0 . Karma and Rappel show how their approach allows 
the effects of kinetics to be eliminated in the standard phase-field model by an appropriate 
choice of the parameters in the model. We also address this issue for the generalized phase 
field model, and find that the effects of kinetics can be eliminated for a significant range of 
the ratio of thermal conductivities as well. 

In the next section we describe the phase field models that we consider. Asymptotic 
analyses of these models are presented in Section 3, followed by discussion and concluding 
sections. 

2 Phase-Field Model 

A common version of the standard phase-field equations for a pure material can be written 
in the form 

T* = KV 2 4> - y g'(4.) - (2) 

c P T t - L v r'{<!,yp t = V • [mVT] , (3) 

where the phase field (f>{x , t) is a nonconserved order parameter that labels the phases, with 
0—1 in the solid phase and 0 = 0 in the liquid phase. Here M > 0 is a mobility coefficient, 
K > 0 is a gradient energy coefficient, W > 0 is a barrier height for the double well 
potential function g(<fi) = 0 2 ( 1 — 0) 2 and p(4>) = 0 2 (3 — 20). The function k{4>) permits 
different thermal conductivities ki = k(0) and ks — &(1) in the bulk liquid and solid phases, 
respectively; specific forms for k{4>) will be discussed below. The function r(0) is assumed 
to have r(0) = 0 and r(l) = 1. Examples include the cases r(0) = p(0) and r(0) = 0. 
Following Almgren we refer to the system with r(0) = p(0) and r(0) — 0 as the gradient 
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equations and the nongradient equations, respectively; the reason for this terminology lies 
in how the equations are derived (see [2] and below). 

An isothermal version of this model has an associated bulk free energy density of the 
form 

m, T) = y g{4>) + — T ~ Tm W )- (4) 

l 1 m 

In this case the phase-field equation (2) admits a one-dimensional traveling wave of the form 

<t>( x » t) = ^ j 1 - tanh | , (5) 


where the interface velocity is given by V = 6M LyUT^ — T)/Tm ■ The interface width L 
surface energy 7, and kinetic coefficient 11 that are associated with this solution are given by 


e = y/K/w, 


1 


vkw 


ML v l 

ft rp 

lM 


( 6 ) 


where the latter expression follows from the definition V = /i(Tm — T). Karma and Rappel 
[1] provide an alternate interpretation of the kinetic effect through their thin interface limit, 
which is also discussed below; this limit is intrinsically non-isothermal, so that the above 
interpretation does not apply. 

The governing equations (2) and (3) can be derived from an entropy functional following 
the formalism of irreversible thermodynamics [17, 18]. A generalized version of these equa- 
tions can be derived by using both an internal energy functional S and an entropy functional 
S , which we write in the form [19-24] 


E = J t T) + ^\V^f}dV, 

( 7 ) 


( 8 ) 


where K E and K$ are gradient coefficients that will be assumed to be constant, and e 
and s are the internal energy and entropy densities, respectively. By starting with a free 
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energy density /(T, <f>), it is thermodynamically consistent to take s(T,<f>) — —dfjdT and 
e (T,<l>) — f{T,(j>) + Ts(T,4>). 

An associated free energy functional T can be defined for the system, which is given by 

r=J v {/(<», T) + (Ke ‘J K ' - |V^| 2 | iV. (9) 

Anisotropic versions of the model can be obtained by modification of the gradient terms 
|V0| 2 , as described in more detail below. For simplicity we will assume that the densities 
and specific heats are uniform throughout the system. 

The derivation of the governing equations may be based on conservation laws for energy 
and entropy written in the form 


d£ f =• ^ . . 

— + / Je -ndA — 0, 
dt Jsv 


( 10 ) 


// C p /* 

— + / J s - ndA = / spdV (11) 

dt Jsv Jv 

where Je and Js represent the flux of energy and entropy through the boundary 8V. which 
has outward unit normal n, and sp represents the entropy production term. Governing 
equations that guarantee positive entropy production can be obtained by the choice 


Je — —k,WT — K E <l> t V(l>, (12) 

Js = -^T + K s 4>tV4>, (13) 

where the terms proportional to the gradient coefficients represent non-classical fluxes asso- 
ciated with the interfacial region [18]. The entropy production is then given by 


Ts P - —|VT| 2 + <j> t [(K e + TK S )V 2 (j> - U] , 


(14) 


and the corresponding phase-field equation [cf. Eq. (17)] is chosen so that the latter expres- 
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sion takes the form k\VT\ 2 /T + <j) 2 /M, with M > 0. 

In this paper we consider a free energy of the form 


M T) = {We + J Ws) g(<fr) + Lv(T Tm) pW - cp(Tln[T/Tu\ - T). (15) 

Z ± M 

The corresponding internal entropy and energy densities may be found, using the relations 
s = —df/dT and e = / + Ts, to be 

Wr W<? Lu 

e(</>,T) = cpT+-£g^)-L v pW, = c P ln[T/T m ]-- 1 f-g(4’)--^p(d,). (16) 

/ / 1 M 


Here We and W$ are constant double well heights for internal energy and entropy, respec- 
tively. The phase-field models considered in Refs. [17] and [18] correspond to a specific 
entropic version formulation obtained from our model when Ke — We = 0. 

The generalized phase-field equations corresponding to this model take the form 


i-A, = (K e + TKs)V v ~^ E + TWs) - Lv(T ~ Tm) 


1 

M 


2 


-m - 


T, 


M 


P'(<t > ), (17) 


c P T t - L v p' {(j>)4>t + = V • [k{4>)VT] + K E (j) t V 2 (j). 


(18) 


This generalized model includes two additional temperature-dependent terms (T—Tm)Ks'^7 2 4 > 
and (T—TM)Wsg'(<j >)/ 2 in the phase-field equation, as well as two additional terms W E (j>t g'{4>)/ 2 
and K E <t>t V 2 (£ in the energy equation, when compared to the standard version (2) and (3). 

If these terms are omitted, these equations are similar to the standard version, but with 
p'(<j>) appearing in place of r'{<f>) in the energy equation. We will refer to equations (17) and 
(18) with p((p) replaced by r{(fi) in (18) as the generalized phase- field equations. The forms of 
the generalized phase-field equations corresponding to the particular cases r((f)) — p((j>) and 
r{4 > ) — (j) bo termed the generalized gradient equations and the generalized nongradient 
equations , respectively. 

The question of whether K E , Kg, W E and Wg represent measurable physical quantities 
is an important issue that concerns the physical validity of Eqs. (17) and (18). Here we 
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adopt a pragmatic approach, and view the formal derivation of these equations as a means 
of suggesting modified forms of the governing equations that can be advantageous from the 
viewpoint of asymptotics and computations. For example, Karma and Rappel [1] note that 
there are computational advantages to using nongradient governing equations that do not 
result from a formal thermodynamic derivation; they suggest the use of the function r(< jf>) = © 
rather than r(<f) — p(4>) in the latent heat term of the energy equation. In a similar vein, 
Almgren considers the question of choosing the specific form of the interpolation functions 
p(<j)), r(0), and k((j>) to eliminate unwanted terms in the thin interface asymptotics that arise 
when k,L ^ ks- Here we examine the effect on the thin interface asymptotics of introducing 
the terms involving K E , K s , Wp and Ws for fixed forms for the interpolation functions. 

2.1 Dimensionless Equations 

We make the governing equations (17) and (18) dimensionless by using the macroscopic 
length scale R, the thermal time scale R 2 /R, and the temperature scale L v jc P ; here R = 
(ke + kl)/2 is the average of the thermal diffusivities Kp = kp/cp and k,s = ks/c P in the 
liquid and solid, respectively. We retain the original definitions of the interface width l and 
surface energy 7 by using K — (Kp + TmK s ) and W — {Wp + TmWs ) in Eq. (6). The 
resulting dimensionless equations can be written in the form 

e 2 T<j> t = e 2 (l + au)V 2 4> - (19) 

Ut - r\ 4 >)<t>t + - v • [Q(<f>y Vu] + e 2 ^ t V 2 0, (20) 

where we have defined the conventional dimensionless parameters 

= — A = — _ RT m c P 

e R Uc T Mt c ML 2 v 1 j 
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that usually appear in the phase-field equations, as well as the parameters 

K s {L v /cp] „ W s [LvIcp\ r W E K n 

“= K'~"’ 0= W ' i= I7’ V= LjP' 

that appear in the generalized model. For now we will treat a, j3, S, and v as independent 
parameters, and will return to the question of the relationship between these parameters 
below. We have also defined the dimensionless function Q{<j>) = k(<fi)/[Rcp} with Q(0) = 
Ql = kl/k and Q(l) = Qs — «s/«- We note that since Q L and Qs are normalized in terms 
of their mean value, we have Ql + Qs = 2 and hence 0 < Ql < 2 and 0 < Qs < 2. 

The relative sizes of the dimensionless parameters determine what type of sharp interface 
analysis is appropriate in deriving the limiting form of the associated free boundary problem. 
The original approach to the sharp interface analysis, which we will refer to as the “classical 
analysis,” and the thin interface analysis of Karma and Rappel differ particularly in the 
assumptions concerning the scaling of the factor A u in Eq. (19). 

3 Asymptotic Analyses 

The various sharp interface limits use the method of matched asymptotic expansions, whereby 
an outer solution, valid away from the interface, is matched to an inner solution that is valid 
in the interfacial region [6]; for brevity wo only sketch the procedure. The inner expansion 
near a point Xq on the interface is described via a stretched coordinate £, with x — xo + 
where n is the local unit normal to the interface. The low order terms in an inner expansion 
of the solution can be obtained by using the relations 

v 2 = + o(i). 4 = ^{ + 0(1), (23) 

where the dimensionless mean curvature % = RH and the dimensionless interface speed 
v n — VR/R are assumed to be of order unity. In the inner region the temperature and phase 
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field are expanded in powers of e, 


0 = 0(0) + € 0(D +0 {e 2 ), 


(24) 


u = + eu^ + 0(e 2 ), 


(25) 


and the resulting equations are solved term by term in e, with far-held boundary conditions 
that are obtained by matching to the outer solution. The outer solution has 0(f) equal to 
zero or one to all orders in the liquid or solid respectively, and the thermal field has the 
expansion U(x) = U^\x) + eU^(x) + 0(e 2 ). The limiting behavior near a point Xo on the 
interface is given by 


U (£„ + e(n) — U ± ^ (xo) + e 


U±\xq) + C I +< ^( e2 )’ 


dUf ] 


(26) 


where and denote the limits of U as ( — > 0 + and £ — > 0 _, respectively, as so forth. 
The inner solution must match with this behavior as £ — > ±oo [6]. 


3.1 Classical Analysis 

In this setting a standard model given by (19) and (20) with a, ft, 5, and u all set to zero 
is employed. The asymptotic analysis with i « l c assumes that A — £'/(60) is a small 
parameter, which allows one to take u to be of order unity while relegating to 0(e) the term 
Xu. If we write A — Ae, then the leading order phase-field equation is given by 

< - = 0 , (27) 


with the solution 


* <0) (O = 5 


1 — tanh 



(28) 
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The leading order temperature equation is 

= °- ( 29 ) 

Integrating and matching with the outer solution twice in succession gives that — 
U+\x o) = £/i 0) (£o) is constant. 

The first order phase field equation is then 

4 $ - ^/(^ (O) )0 C1) = -('TVn + + Au ( °V (<£ (0) )> (30) 

and a solvability condition for this equation then gives (c.f. Eq. (38) below) 

um = w-JK = -^ H - v ' v - (31) 

which recovers the modified Gibbs-Thomson boundary condition (1), where the kinetic term 
agrees with the form of V 0 — n{L v /cp\ given in Eqs. (6) and (21); thus both the isothermal 
traveling wave solution and the classical asymptotic analysis lead to the same definition of 
H in terms of the phase field parameters. 

The first order temperature equation is then 

u n r'(^°))4 0) = (g(^°))74 1) ) c . (32) 

By integrating this equation over the interval — oo < £ < oo and using the matching con- 
ditions, the leading order heat flux boundary condition is recovered in the dimensionless 
form 

-Vn = QlGl - QsGs , (33) 

where Gi = OU^ (On and Gs = dU^ / dn denote the local temperature gradients at the 
interface. 

Almgren has considered the forms of the terms that occur at the next order in the classical 
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analysis, and has discussed strategies for eliminating anomalous terms at that order so as 
to obtain higher-order accuracy in the Gibbs-Thomson equation [2]. Almgren also discusses 
the thin interface analysis of Karma and Rappel [1], to which we turn next. 

3.2 Thin Interface Analysis 

In the thin interface, or isothermal, analysis, A is assumed to be of order unity, but Aw = 
0(e) is still retained by taking u ^ = 0 so that u = 0(e). We will show that for the 
generalized phase-field equations the quantities a, ft, S and v may be chosen to ensure that 
both the interfacial temperature is continuous and the modified Gibbs-Thomson equation is 
independent of Gl and Gs- 

With the assumption that the parameters A, r, v n , H, a, ft, S and v are all of order unity, 
the leading and first order phase-field equations are then 

4>a ~ = 0, (34) 

4 ) -^ (0) )^ (1) = i? i’ ( 35 ) 

where 

Ri = —au^4^ + ^w^V(^) — ( TV n + R)4 > f > + (36) 

The leading order phase field solution is again given by Eq. (28) . The first two terms of Ri 
given by Eq. (36) can be simplified using Eq. (34) to give 

-aw (1 V[c + ^u (1 V(0 (o) ) = | u (1 V(^ (0) )i (37) 

where we have defined ft — ft — a. The solvability condition for Eq. (35) can then be written 
in the form 

4>f ] Rid( = A f w (1 V(^ (0) )4 dC + ? / u (i y(^> (0) )<^ 0) d( - ~(rv n + H), (38) 

-oo J — oo Z J — oo 0 
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which involves the first order temperature solution We observe that only the combina- 
tion /? — a — j3 enters the solvability condition. 

With = 0 the temperature equation is satisfied identically to leading order. The first 
order temperature equation is 

v n r'(0 (o) )^ o) - ^v n y ' («^ (0) )^ 0) = ^ - WntfUft-, (39) 

which, when simplified by using Eq. (34) as above, may be integrated to yield 




u 


( 1 ) 


V n T 




2 


1 2 


= QlGl — QsGs 


(40) 


where we have defined v = v — 8. Here the constant of integration has been evaluated in 
two ways by taking £ — > ±oo and applying the matching conditions (26). Note that only 
the combination v — 8 = v appears in this expression. The solution can then be expressed 
in the equivalent forms 


u (l) (0 = 

- l°° opsjy {vW <0) ) + ^ [4 0) ] 2 - G l [Q(^) - Ql]} dr, (41) 

= U"’ix„) + CG S 

+ /_£ Q^yj {«»K* (0) ) - 1] + *~y K“T - G S [W <0) ) - <?*]} dr,. (42) 

3.2.1 Continuity of Temperature 

Subtracting the two equivalent forms for m^O) provides an expression for the jump in the 
temperature across the interface, 

VftHi*) - V m (fo) = l °° opijy [4 0) ] 2 - Gl[W w ) - Q,.]} dr, 

+ /_£ - 1] + X [«f ] 2 - GslQ(^°>) - Qs)}d„. (43) 
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For v = 0, Karma and Rappel [1] found that this expression vanishes for either ~ 

<f> or r((f>) = p(4>) if the thermal properties in the liquid and solid phases are the same. 
Almgren [2] considered the more general case of unequal thermal properties, and noted that 
it is convenient to interpolate 1 /Q{4>), since that quantity appears in the above integrand. 
Almgren discusses how to eliminate possible temperature discontinuities by choosing more 
elaborate interpolation functions for the case with V = 0. Here for V / 0 we consider the 
simple choice 


= [1 ~ rW) . . 

W) Qs Ql ' ^ ’ 

where we interpolate using r(<f>). For r (<j>) — <f> the integrals in Eq. (43) can then be evaluated 

directly to yield 


u +\%o) - {s + 




(45) 


and for the choice r((j>) = p{4>) — (f > 2 ( 3 — 2(f) we similarly find 


U^(x o) - U?\x, o) = v n 


v_ ( J_ J_\ 19 / J_ 

24 + + 30 



(46) 


In performing the integration, we note that the terms involving Gl and Gs in the integrands 
of Ecp (43) combine and simplify upon using the heat flux condition (40). Thus the leading 
order temperature field is continuous across the interface if we choose 


v = 12 (Q l - Q s ) = 24 


f KL ~ Ks \ 
\KJL + KsJ 


in the case r{(f) = (j), and choose 


v — 



(47) 


(48) 


in the case r{(f) = p(4>)- Values of v for other choices for r(<f) can be calculated in a similar 
manner to ensure that the interfacial temperature is continuous. Almgren [2] discusses 
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strategies for the standard model to ensure that the usual boundary conditions are satisfied 
with 0(e 2 ) errors as well. 


3.2.2 The Modified Gibbs-Thomson Equation 

We now investigate the solvability condition (38) to determine the form of the modified 
Gibbs-Thomson equation in the thin interface limit. Evaluating the integrals in Eq. (38) 
and using the relation 


G c -G l = 


J_ 1 

Ql Qs 


-G 


2 i_ 

Ql Qs 


(49) 


where G = ( QlGl + QsGs)/ 2 is the (conductivity-weighted) average temperature gradient, 
gives that the solvability equation can be -written in the form 


-(c4 1, + c <1) ) = »„((— + -E) 

2^ + t n \ \ 560 240/ 


_L j_ 

Ql Qs 


5A fiv \ 




0 T 


_1 1_ 

Ql Qs 


+ 24 


2_ _i_ 

Ql Qs 


12 1120 y 

) ' 


1 1 
Q~l + Q~s 


n 

6 ’ 


for r((/>) = <f>, and 

1 ([/'■> + e (I) ) =vJ(— + -G 

2 V + > ” 1 V 144 144. 


1 1 

Ql Qs 


+ 


19A _ pv A 


1 1 

~60 1120 J [Ch + Q~ s 


T 

6 


(50) 


T 

6 


G 


r i9A 

' 1 1 ' 

| 30 

Ql Qs. 


+ 24 


1 1 
Ql Qs 


n 

6 ’ 


for r(4>) = p{4>). 

The term proportional to G is eliminated if we choose 


(51) 


P = 10A (Q l - Q s ) = 20A 


(52) 


15 



in the ease r(<f>) = <j>, and 


fi = f A (Ql ~ Qs ) = j A 


f ~ K S 
\Kl + Ks 


) 


(53) 


in the case r(<^>) = p(<f>). 

Inserting the expressions (52) or (53) for ft and (47) or (48) for v into the solvability 
condition then gives 


6A U w = -v n 



CjA • 

QlQs - 


1 — C 2 (Ql — 



-n 


(54) 


where C\ and c 2 are purely numerical constants whose values are given by Ci = 5 and 
c -2 — 269/700 for r(tf>) = <j>, and ci = 19/5 and c 2 = 173/525 for r(4>) — p(4>). This expression 
can be written in terms of dimensional variables as 


u = eU {1) = -V/Vq - 4 H, 


(55) 


where Vo = p[Ly/cp\ and 


1 _ 

fj, 6£ML V 


ci£k[Lv / cp\ 

Qk l ks 


1 - 4c 2 


( kl- ks 
\k l + k s 



(56) 


Following Karma and Rappel [1], we note that the effects of kinetics can be eliminated 
(l//j, = 0) by the choice 


1 _ a^K^/cp] 1 _ 4c / *l ~ K s \ 1 2 
M kl^sTm 2 \kl + k,s) 


(57) 


In order that the phase-field equation is well posed we require that M > 0 and so (57) may 
be satisfied providing that the ratio Ks/kl is not too extreme; specifically 


— 1 «£ + 1 
2-y^+l k l 2^-1’ 
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which evaluates to 0.107 < ks/k>l < 9.34 for r{<j>) = cf) and 0.0689 < ks/ k l < 14.5 for 
r(4>) = p(<j>). Both of these criteria are met by most semiconductors and metals as well as 
succinonitrile. 


4 Discussion 

In this section we address a number of related issues that concern the asymptotic analysis. 
Wc first, discuss the elimination of anomalous terms by using the generalized phase-field 
model, followed by an indication of how surface tension anisotropy can be included in the 
model. The magnitudes of the anomalous terms in the thin interface analysis of the standard 
phase-field model are then estimated. We conclude the section by comparing the assumed 
scalings for the classical and thin interface asymptotics. 

4.1 Elimination of Anomalous Terms 

In the isothermal analysis, we find that of the four parameters a, /3, 5, and v that occur in the 
generalized phase field equations, only two parameters, ft — ft — a and v — v — 6, appear as 
degrees of freedom that are available in the first order analysis to ensure that the interfacial 
temperature is continuous and that the anomalous term that is proportional to the 
temperature gradient in the Gibbs-Thomson equation vanishes. If, instead of assuming the 
independence of the dimensionless parameters a, ft, 6, and v we consider their definitions 
in terms of the parameters K E , Ks, We and Ws that appear in the internal energy and 
entropy functions, we find that (3 and V are actually not independent, but are related, since 
Eq. (22) implies that 

P - AP - --p—(W E K s - W s K e ). (59) 

6 bcpj^ 

Thus, if the parameters in the generalized phase-field equations have values that are derived 
from an underlying internal energy and entropy functions in a thermodynamically-eonsistent 
manner, only a single degree of freedom, namely W E Kg — WsK E , is available to satisfy the 
two conditions. However, we note that for the generalized gradient equations , these conditions 
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given by (48) and (53) are linearly independent and reduce to 


P = \v = \v\ 


( 60 ) 


where 




76 (k l - k s 


5 \k, l + k s 


(61) 


Hence, in the fully thermodynamically consistent treatment, represented by the generalized 
gradient equations, we can choose the four parameters a, ft, 6, v so that both the tempera- 
ture at the interface is continuous and the Gibbs-Thomson equation is independent of the 
interfacial temperature gradients. In addition, for a wide range of values of ks/kl M may 
be chosen using equation (57) so that kinetic effects are not present, i.e., the dimensional 
boundary conditions on the value of the interfacial temperature are 


T s = T l = T m - 

Cp 


(62) 


where Ts and Tl are the interfacial temperatures in the solid and liquid respectively. We 
observe that the choice of a, /3, S, v consistent with (60) is not unique. For example, we 
may choose a — 5 — 0, f3/\ — v = V*, which corresponds to allowing the bulk free energy 
double well height W = We + TWs to depend on T and the gradient energy coefficient 
Kf — Ke + TKs to be independent of T. Similarly, the choice /? = v = 0, a/A = S = — v* , 
corresponds to the opposite situation where W is independent of T and K F is dependent on 
T. This latter choice may be advantageous in numerical computation using the generalized 
gradient equations since it eliminates the 0tV 2 ^> term in the Eq.(20). In the nongradient 
case it is still possible to eliminate the anomalous terms using the choices for ft and v given 
by equations (53) and (48), respectively; however, in this case the condition between them 
given by equation (59) will not be satisfied. 

The isothermal asymptotics applied to either of the generalized phase-field equations 
gives interfacial conditions that depend on the forms adopted for the functions r(<j)), Q(4>) 
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and p(4>). In particular, the results given above for the generalized gradient equations are 
for the specific choice p(0) = 0 2 ( 3 — 20) and the form of Q((f>) given by equation (44). We 
conjecture that for a more general form of p(<f>) the elimination of the anomalous terms 
will still require the relationship (3 = XV to hold but will alter the expression for v*. In 
contrast, the classical asymptotics gives interfacial conditions on the temperature that are 
independent of the forms of r(0), Q(0) and p(4>). 

In the thin interface limit, the sensitivity of the limiting interfacial boundary conditions 
to the details of the model within the interface can be attributed to the relatively large 
interfacial thickness allowed by this limit as compared to the classical sharp interface limit. 
A similar sensitivity was also observed in a different thin interface analysis of solute trapping 
for a phase field model of a binary alloy [28] . In this situation the dependence on the details 
of the phase-field model in the interfacial region was apparent in an asymptotic regime 
where the interface thickness was also comparatively large, in fact comparable to the solute 
boundary layer thickness. 

We note that for both the classical analysis and the thin interface analysis the temperature 
field near the interface is spatially uniform to leading order in e; in the thin interface analysis 
T = Tm to leading order, and in the classical analysis there are corrections to T M of order 
unity due to capillary and kinetic terms. The first order correction to the temperature field 
is non-constant in each analysis, with an 0(e) spatial variation in temperature through the 
interface. 

4.2 Anisotropic Surface Energy 

The phase field equations (19) and (20) are based on isotropic forms of the gradient energy 
and gradient entropy represented by the |V0| 2 terms in the functionals (7) and (8). A 
model with anisotropy can be obtained by replacing |V0| 2 in these expressions by the term 
[r(V0)] 2 , where T(p) is a homogeneous function of degree one [8,25-27]. The introduction 
of the function T(V0) provides a phase-field model of interfaces whose anisotropic surface 
free energy is proportional to T(n) in the sharp interface limit, where fi is the outward unit 
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normal to the interface. For example, the function 


r (V0) = |V0| {1 + 74 cos (4arctan {(j) y l<t>x])} , (63) 

models a two-dimensional interface with a four-fold surface energy proportional to T(n) = 
(1 + 74 cos 40), where 9 is the orientation of n — (cos 6, sin 9). The resulting anisotropic forms 
of the generalized phase-field equations are obtained by the substitution 

V 2 0 -+ V • [r^ (64) 

in equations (19) and (20), where £(V0) is the Cahn-Hoffman ^-vector [30] with components 
£j{p) — dT (p)/dpj. The sharp interface analysis of the anisotropic equations parallels the 
treatment in Ref. [25]; for brevity, we omit the details. The final result is that the mean 
curvature term H in the Gibbs-Thomson equation (55) is replaced by the weighted mean 
curvature term Vs*£ involving the surface divergence of the ^-vector, which is the appropriate 
generalization of the Gibbs-Thomson equation for the anisotropic case [31,32]. 

4.3 Magnitudes of Anomalous Terms in the Thin Interface Limit 

In Section 3.1 we sketch the classical analysis for the standard model. It is straightforward 
to perform the classical analysis for the generalized model as well. One finds that the leading 
order temperature is continuous at the interface and satisfies a non-linear version of the usual 
(modified) Gibbs-Thomson equation. The conventional leading-order balance of heat flux at 
the interface is also obtained, but with a discontinuity in the first-order temperature field 
across the interface. It is this observation that suggests the use of the generalized phase-field 
model in the thin interface limit in this paper to eliminate the temperature discontinuity 
that arises for unequal thermal conductivities. 

The thin interface limit of the conventional phase-field model also leads to a jump in 
temperature at the interface given by Eq. (45) or (46) with i> — 0. Considering as an 
example the case — <j>, the dimensional form of the leading order temperature jump can 
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be written as 

Tl ~ Ts = lLv/Cp] (* ~ fe) yr L . ■ (65) 

which shows that the temperature jump vanishes if the thermal properties of liquid and 
solid are equal. The jump in temperature is a non-equilibrium effect that resembles in some 
ways the “solute trapping” effect that is observed for rapid solidification of a binary alloy 
[33]. Solute trapping involves non-equilibrium partitioning in the solute at the solid-liquid 
interface due to a jump in chemical potential across the interface. The jump in potential 
becomes significant at velocities V ~ Vp, where the characteristic velocity Vp is given by 
the ratio of an interfacial solute diffusivity Dj to the interface width l. In a similar way, 
the temperature jump in (65) is proportional to V/Vp, where Vp = is a characteristic 
thermal velocity. In phase-field models of solute trapping, Vp is found to be sensitive to 
the assumed form of the variation of the solute diffusivity through the interface [28], just 
as the numerical factors in Eqs. (45) and (46) depend on the assumed form for r ((/)). The 
sensitivity of the limiting interfacial boundary conditions to the details of the model within 
the interface in the thin interface limit can be attributed to the relatively large interfacial 
thickness allowed by this limit, as compared to the classical sharp interface limit. For 
solute trapping the appropriate limit also involves an asymptotic regime where the interface 
thickness is relatively large, in fact comparable to the solute boundary layer thickness. For 
solute trapping in metallic alloys a typical value of Vp is on the order of 1 m/s; since the 
thermal diffusivities are usually orders of magnitude larger than typical solute diffusivities, 
the temperature jump is a much smaller effect. Using values for lead as an example with 
ks/k L w 2, a nominal interface thickness of 10 7 cm leads to a temperature jump of about 
10~ 2 K at a velocity of 1 m/s. 

The thin interface limit also results in an anomalous term in the modified Gibbs-Thomson 
equation, which for the case r(4>) = <j) can be written in the dimensional form 
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where the last term in this expression depends on the conductivity-weighted temperature 
gradient at the interface. This term is a very small correction for metallic systems under 
most conditions; for example, using values for lead with a mean temperature gradient of 100 
K/cm and £ = 10 -7 cm gives a correction of roughly 10~ 5 K. 


4.4 Scalings of Surface Energy and Kinetics 

We now return to the question of the assumed scaling of capillary and kinetics effects in the 
classical and thin interface asymptotics for the standard equations. We couch the discussion 
in terms of the dimensionless parameters e = £/R, A = £/(6£ c ), and r, where we recall that 
£ is the interface width, R is a macroscopic length scale, £ c = Tm') /[L y / c p] is the capillary 
length, and r is a dimensionless form of the mobility coefficient M [cf. Eq. (21)] that is 
assumed to be of order unity in both analyses. The two asymptotic limits correspond to 
e — ¥ 0 with either A — 0(e) (the classical analysis) or A = 0(1) (the thin interface analysis). 
The gradient energy coefficient K, double well height W , and mobility coefficient M can be 
written in the form 




W = 6 


— ) 

cpT M ) 


1 

A’ 


/ kTmCp \ 

\t L yR 2 ) e 2 


(67) 


In the classical analysis the surface energy and kinetic coefficients have the form 


( L 2 R \ e 1 _ ( tLvR \ e 

\c p T m ) A’ fi V 6 cpR ) A 


(68) 


In the thin interface analysis, 7 also has this form, but 1/ \i has an additional term as given 
by Eq. (56). In the classical analysis with A = O(e), 7 and 1/p, scale as 0(1) quantities with 
respect to e, the gradient energy coefficient is 0(e), the double well height is 0(l/e), and 
the mobility is 0(l/e). In the thin interface limit with A = 0(1), 7 tends to zero as O(e), 
the gradient energy coefficient tends to zero as 0(e 2 ), and the double well height scales as 
an 0(1) quantity. The interpretation of the kinetic coefficient in this limit is more involved, 
since one can arrange that l//r = 0 by an appropriate choice of M as given by Eq. (57). 
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With this choice, M scales as 0(l/c 2 ) which is consistent with the scaling given by Eq. (67) 
with A — 0(1). This shows that, in the thin interface limit, the effects of capillarity and 
kinetics are assumed to be small corrections in the analysis, as noted in the introduction. 

In the context of numerical calculations using a phase-field formulation, the two asymp- 
totic analyses suggest different interpretations of a given computation, and different strategies 
for the appropriate selection of parameters for a sequence of calculations that are intended 
to reproduce a sharp-interface model. With the classical analysis, to converge to a sharp 
interface result one should in principle perform a sequence of calculations with £ tending to 
zero for fixed values of 7 and 1 / p, so that e and A both tend to zero. By contrast, for the 
thin interface analysis, the surface energy and kinetic effects vanish as £ tends to zero, and 
in principle a single calculation with a fixed (small) value of e is performed, with appropriate 
values prescribed for K, W, and M in order to give the desired values of 7 and 1 / /i consistent 
with the chosen value of e. Performing a convergence study with e tending to zero in the thin 
interface limit requires constant reduction of A to preserve given values of surface energy and 
kinetics, causing the procedure to revert to the classical analysis in the limit of small e. On 
the other hand, lower bounds on the feasible size of e are set by computational constraints 
in practise, so convergence studies of this nature have limited applicability. 

5 Conclusions 

The thin interface asymptotics of Karma and Rappel applied to the standard phase-field 
equations for a pure material with unequal thermal conductivities generate anomalous terms 
involving a temperature jump across the interface and a contribution from the mean temper- 
ature gradient in the modified Gibbs-Thomson equation [1,2]. The temperature discontinuity 
can be interpreted as a thermal trapping effect analogous to the solute trapping that occurs 
during the rapid solidification of a binary alloy. For metallic systems under typical growth 
conditions both anomalous terms are very small, and can be eliminated by employing a 
generalized phase-field model that is based on entropy and internal energy gradient energy 
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functionals that include double wells in both the internal energy and entropy densities. The 
flexibility of the thin interface analysis in eliminating kinetic effects carries over to the case of 
unequal thermal conductivities as well, as also noted by Almgren in his isothermal analysis 
of this limit. In our approach kinetic effects may be eliminated for conductivity ratios of up 
to about 10:1. The use of thin interface asymptotics to eliminate kinetic effects in binary 
alloy systems where diffusivities may differ by several orders of magnitude is an area of active 
research, where the use of generalized energy-entropy models may play a useful role. 
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